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We study the optomechanical design introduced by M. Li et al. [Phys. Rev. Lett. 103, 223901 
(2009)], which is very effective for investigations of the effects of reactive coupling. We show the 
normal mode splitting that is due solely to reactive coupling rather than due to dispersive coupling. 
We suggest feeding the waveguide with a pump field along with a probe field and scanning the 
output probe for evidence of reactive-coupling-induced normal mode splitting. 
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I. INTRODUCTION 

In a recent paper Li et al. [l| presented a new de- 
sign for an optomechanical system that consists of a mi- 
crodisk resonator coupled to a waveguide. This design 
has several attractive features. Besides its universality, it 
enables one to study the reactive effects [l|, in optome- 
chanical coupling. The origin of the reactive coupling is 
well explained in Ref. Q. Its origin lies in the mechani- 
cal motion dependence of the extrinsic losses of the disk 
resonator. Further phase-dependent gradient forces lead 
to reactive coupling. Li et al. have also argued that this 
design is more effective in achieving cooling of the sys- 
tem to its ground state. While cooling is desirable for 
studying quantum effects at the macroscopic scale 
we examine other possibilities, which do not depend on 
the cooling of the system, to investigate the effects arising 
from strong reactive coupling. Since optomechanical cou- 
pling effects are intrinsically nonlinear, we examine the 
nonlinear response of the microdisk resonator to pump 
probe fields. We report reactive-coupling-induced nor- 
mal mode splitting. Note that in previous works [lfj| - fl~3| 
on normal mode splitting in optomechanical devices, only 
dispersive coupling was used. In this paper, we report on 
normal mode splitting due to reactive effects. 

The paper is organized as follows. In Sec. II, the phys- 
ical system is introduced and the time evolutions of the 
expectation values of the system operators are given and 
solved. In Sec. Ill, the expectation value of the output 
fields is calculated, and the nonlinear susceptibilities for 
Stokes and anti-Stokes processes are obtained. In Sec. 
IV, we discuss normal mode splitting in output fields 
with or without reactive coupling. We find that there is 
no normal mode splitting in output fields in the absence 
of reactive coupling. However, normal mode splitting oc- 
curs in output fields in the presence of reactive coupling. 



II. MODEL 

We consider the system shown in Fig. [T] in which 
a microdisk cavity is coupled to a freestanding waveg- 
uide. A strong pump field with frequency u>i and a weak 



Stokes field with frequency lu s enter the system through 
the waveguide. The waveguide will move along the y di- 
rection under the action of the optical force exerted by 
the photons from the cavity. Further, considering the 
dispersive coupling and reactive coupling between the 
waveguide and the cavity, displacement q of the waveg- 
uide from its equilibrium position will change the reso- 
nant frequency of the cavity field and the cavity decay 
rate, represented by uj c (q) and K e (q), respectively. 




FIG. 1: (Color online) Sketch of the studied system (from 
Ref.l). The microdisk cavity is driven by a pump field and a 
Stokes field. The nonlinearity of the interaction also generates 
anti-Stokes field. 

In a rotating frame at pump frequency w/, the Hamil- 
tonian of the system is given by [l| 
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The first term is the energy of the cavity field, whose 
annihilation (creation) operators are denoted c(c T ). The 
second and third terms are the energy of the waveguide 
with mass m, frequency u) m , and momentum operator 
p. The fourth term gives the interactions between the 
waveguide and the incident fields (the pump field and 
the Stokes field), L is the length of the waveguide, c is 
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the speed of light in vacuum n g is the group index of the 



waveguide optical mode [14| . £; and |e s | are the ampli- 
tudes of the pump field and the Stokes field, respectively, 
and they are related to their corresponding power pi and 
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describe the coupling of the cavity field to the pump field 
and the Stokes field, respectively. And 5 = u s — cj ; is the 
detuning between the Stokes field and the pump field. We 
would study the physical effects by scanning the Stokes 
laser. 

For a small displacement q, w c (<?) and K e (q) can be 
expanded to the first order of q, 



uj c (q) « u c + qx, 

K e{q) ~ K e + q^omi 
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thus the quantities \ an d K om describe the cavity- 
waveguide dispersive and reactive coupling strength, re- 
spectively. Further, note that the photons in the cavity 
can leak out of the cavity by an intrinsic damping rate 
Ki of the cavity and by a rate of K e {q) due to the reac- 
tive coupling between the waveguide and the cavity. In 
addition, the velocity of the waveguide is damped at a 
rate of 7 m . Applying the Heisenberg equation of motion 
and adding the damping terms, the time evolutions of 
the expectation values ((<?), (p), and (c)) for the system 
can be expressed as 
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Eq. ((4]) into Eq. ([3]), ignoring those terms containing the 
small quantities e 2 s , e* 2 , |e s | 2 , and equating coefficients 
of terms with the same frequency, respectively, we obtain 
the following results 
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where 



A = uj c - ui + XQo, 



d(S) = V*(BE + FJ) + BF*J*, 
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and A = \/k{1 -I — —Qo), B = k + qoK om — £(A + 5), 
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approach used in this paper is similar to our earlier work 
13] which dealt with optomechanical systems with dis- 
persive coupling only. 
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III. OUTPUT FIELDS 
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where we have used the mean field assumption (qc) = 
(q)(c), expanded n e {q) to the first order of q, and assumed 
K e rs Ki ~ k/2, where k is the half-linewidth of the cavity 
field. It should be noted that the steady-state solution of 
Eq. ([3]) contains an infinite number of frequencies. Since 
the Stokes field e s is much weaker than the pump field 
the steady-state solution of Eq. ((3]) can be simplified 
to first order in e s only. We find that in the limit t — > oo, 
each (q),(p), and (c) has the form 
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where s stands for any of the three quantities q, p, and c. 
Thus the expectation values ((<?), (p), and (c)) oscillate 
at three frequencies (uji, uj s , and 2u>i — uj s ). Substituting 



To investigate the normal mode splitting of the output 
fields, we need to calculate their expectation value. It 
can be obtained by using the input-output relation [TBj 
(c out ) = y/2K e (q)(c). If we write (c out ) as 
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where c/ is the response at the pump frequency u>i, c s is 
the response at the Stokes frequency w s , and c as is the 
field generated at the new anti-Stokes frequency 2u>i — u> s . 
Then we have 



ci = \/k(1 + ^-q )c , 
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FIG. 2: (Color online) The real roots of d(S) in the domain 
Re(<5) > as a function of the pump power pi for K om — 
(dotted curve) and n om = — 2ir x 26.6 MHz/nm (solid curve). 
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FIG. 3: (Color online) Imaginary parts of the roots of d(S) as 
a function of the pump power pi for K om = (dotted curve) 
and K om = — 27r x 26.6 MHz/nm (solid curve). 



Furthermore, whether there is normal mode splitting in 
the output fields is determined by the roots of the de- 
nominator d(S) of c s . Here we examine the roots of d(5) 
given by Eq. ([7]) numerically. 

The response of the system is expected to be especially 
significant if we choose lu s corresponding to a sideband 
u s = uji ± u m or uj s — loi ± A, so we consider the case 
A = Lu m . The other parameters are chosen from a recent 
experiment focusing on the effect of the reactive force 
on the waveguide pj: the wavelength of the laser A = 
2ttc/o;; = 1564.25 nm, x = 2tt x 2 MHz/nm, m = 2 
pg (density of the silicon waveguide, 2.33 g/cm 3 ; length, 
10 /im; width, 300 nm; height, 300 nm), n = 0.2io m , 
u) m = 2tt x 25.45 MHz, and the mechanical quality factor 
Q = uJm/lm = 5000. In the following, we work in the 
stable regime of the system. 

Figure [2] shows the variation of the real parts of the 
roots of d(S) in the domain Re(<5) > with increasing 
pump power for no reactive coupling, n om = 0, and for 
Hom = — 27r x 26.6 MHz/nm. For n om = 0, the interac- 
tion of the waveguide with the cavity is purely dispersive; 
the cavity decay rate does not depend on the displace- 
ment of the waveguide. In this case, the real parts of the 
roots of d(6) always have two equal values with increas- 
ing pump power. Thus there is no splitting because the 
dispersive coupling is not strong enough. However, for 



Kom = —2tt x 26.6 MHz/nm, the system has both dis- 
persive and reactive couplings, the cavity decay rate de- 
pends on the displacement of the waveguide, and the real 
parts of the roots of d(S) will change from two equal val- 
ues to two different values with increasing pump power. 
And the difference between two real parts of the roots 
of d(S) in the domain Re((5) > is increased with in- 
creasing pump power. Therefore, the reactive coupling 
between the waveguide and the cavity can result in nor- 
mal mode splitting of the output fields, and the peak 
separation becomes larger with increasing pump power. 
Figure [3] shows the variation of the imaginary parts of 
the roots of d(5) with increasing pump power for zero 
reactive coupling n om = and nonzero reactive coupling 
Kom = — 27r x 26.6 MHz/nm. For n om = 0, the imaginary 
parts of the roots of d{$) do not change with increasing 
pump power. However, for K om — — 27r x 26.6 MHz, the 
imaginary parts of the roots of d(S) change with increas- 
ing pump power. We thus conclude that for the present 
microdisk resonator coupled to a waveguide the normal 
mode splitting is solely due to the reactive coupling. 




FIG. 4: (Color online) The lower two curves show the normal- 
ized quadrature v 3 as a function of the normalized detuning 
between the Stokes field and the pump field, 5/u) m for n m = 
(dotted curve) and n om — —2tt x 26.6 MHz/nm (solid curve) 
for pump power pi — 20 //W. The upper two curves give the 
normalized quadrature w s +1.5 for pump power pi — 200 /jW. 



IV. NORMAL MODE SPLITTING IN OUTPUT 
FIELDS 

We now discuss how the output fields depend on the 
behavior of the roots of d(5). For convenience, we nor- 
malize all quantities to the input Stokes power p s . As- 
suming that e s is real, we express the output power at the 
Stokes frequency uj s in terms of the input Stokes power 
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Further, we introduce the two quadratures of the Stokes 



component of the output fields by v s 



and 
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FIG. 5: (Color online) The lower two curves show the normal- 
ized output power G s as a function of the normalized detuning 
between the Stokes field and the pump field, S/uj m for n om ~ 
(dotted curve) and n om — —2tt x 26.6 MHz/nm (solid curve) 
for pump power pi = 20 /iW. The upper two curves give the 
normalized output power G s + 1.5 for pump power pi = 200 
fiW. 



active coupling (K om —0) and with it (K om = — 2ir x 26.6 
MHz/nm), for different pump powers in Figs. HHH For 
K O m=0, it is found that v s has a Lorentzian lineshape 
corresponding to the absorptive behavior. Note that v s 
and G s exhibit no splitting when n om =0. However, for 
Hom = — 27r x 26.6 MHz/nm, it is clearly seen that normal 
mode splitting appears in v s and G s . Therefore reactive 
coupling can lead to the appearance of normal mode split- 
ting in the output Stokes field. And the peak separation 
increases with increasing pump power [16]. The dip at 
the line center exhibits power broadening. We also find 
that the Stokes field can be amplified by the stimulated 
process. Obviously the maximum gain G s for the Stokes 
field depends on the system parameters. For a pump 
power pi = 200 /^W, the maximum gain for the Stokes 
field is about 1.3. 

Note that the nonlinear nature of the reactive cou- 
pling generates anti-Stokes radiation. In a similar way, 
we define a normalized output power at the anti-Stokes 
frequency 2w; — u> s as 
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FIG. 6: (Color online) The lower two curves show the nor- 
malized output power G as as a function of the normalized 
detuning between the Stokes field and the pump field, 5/cu m 
for K om — (dotted curve) and K om = — 2-7T x 26.6 MHz/nm 
(solid curve) for pump power pi — 20 fiW. The upper two 
curves give the normalized output power G a s +0.15 for pump 
power pi = 200 fiW. 



One can measure either the quadratures 
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of the output by homodyne techniques or the intensity 
of the output. For brevity, we only show v s and G s as a 
function of the normalized detuning between the Stokes 
field and the pump field 5 / uj m for this model, without re- 
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The plots of Gas versus the normalized detuning be- 
tween the Stokes field and the pump field S/uj m for this 
model, without reactive coupling (k oto =0) and with it 
(Kom — — 27T x 26.6 MHz/nm), for different pump pow- 
ers are presented in Fig. [51 We can see that G as ps 
for K om —0. The reason is that the dispersive coupling 
constant x is too small. However, for K om = —2ir x 26.6 
MHz/nm, G as is not equal to zero. This shows that the 
optomechanical system can generate an anti-Stokes field 
with frequency (2u>i — u> s ) due to the reactive coupling. 
For pump power pi = 200 ^W, the maximum gain de- 
fined with reference to the input Stokes power for the 
anti-Stokes field is about 0.1. 



V. CONCLUSIONS 

In conclusion, we have observed normal mode split- 
ting of output fields due to reactive coupling between the 
waveguide and the cavity. Meanwhile, the separation of 
the peaks increases for larger pump powers. Further, the 
reactive coupling can also cause four- wave mixing, which 
creates an anti-Stokes component generated by the op- 
tomechanical system. 
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